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The enhancement of the At, and A c decay rates due to four-quark operators is calculated. Hard 
gluon exchange modifies the wave function for the hu pair in Ab, |\t(Q)|i tt , and the wave function 
for the cd pair in A c , |*(0)|e d . The modified wave function is found to account for the discrepancy 
of 0.20 ps -1 among the lifetimes of A(, and B° and for the difference of 2.6 ps -1 in the case of A c 
and D°. 

Heavy flavour hadrons, H, contain a heavy quark, Q (b and c), and a light cloud comprising of light quarks (anti- 
quarks) and gluons. When the mass of the heavy quark goes to infinity, the picture of the heavy hadron decays are 
so simplified that the light cloud has no role to play. Theoretical study of the heavy hadrons is facilitated by the 
expansion of the hadronic matrix clement, based on the operator product expansion of QCD, in inverse powers of 
the heavy quark mass, m jlj. The leading order of the expansion, corresponding to the asymptotic limit of the heavy 
quark mass, describes the decay rate of a heavy flavour hadron as if it is of a free heavy quark decay. This implies 
the same lifetime for all heavy hadrons of a given heavy flavour quantum number. The next-to-leading order terms, 
appearing at 1/m 2 , describe the motion of the heavy quark inside the hadron and the chromomagnetic interaction, 
, distinguish the lifetimes of the mesons on one hand and the baryons on the other (with exception of ^q). The term 
corresponding to the third order in 1 /m of the expansion is the matrix element of four-quark operators which contains 
effects such as W-exchange, weak annihilation and Pauli interference, coming from the spectator (light) quarks. At 
■ this order, the lifetimes of the various mesons differ among themselves and so is for baryons. The spectator effects 
^ ' which appear through four-quark operators, are supposed, yet poorly-understood, to explain the intricacies of the 
C^ - lifetime differences and to fix the lifetime hierarchy of the hadrons. 

According to the idea of heavy quark expansion, all hadrons of a given heavy flavour are expected to have nearly 
the same lifetime. The theoretical prediction of the ratio of lifetimes of Ah and B°, obtained to two orders in 1/m, is 
0.9. But this is found much higher than the observed value of the ratio: 0.78 for r(Ab) = 1.20 ps and for t(B°) = 1.58 
, ps. On the other hand, the corresponding charm sector shows an intrinsically different picture due to the mass of the 
Q\ ■ charm quark, which is not so asymptotically large as the b quark mass. In the charm case, the dominant effects come 
from the four-quark operators rather than the kinetic and chromomagnetic operators. Experimentally, the ratio of 
the lifetimes of t(A c ) and t(D°) is 0.496 for r(A c ) = 0.206 ps and t(D°) = 0.415 ps. Therefore, for both the cases, 
the explanation for the suspected small lifetime, and an equivalently enhanced decay rate, of Aq should come from 
Qh, the third order term of the heavy quark expansion in 1/m where the matrix element involves four-quark operators. 
Given the present experimental decay rates for beauty hadrons: r(Af,) = 0.83 ± 0.02 ps^ 1 and T(B°) — 0.63 ± 0.05 
ps -1 , the needed enhancement is 0.2 ps -1 , whereas 2.6 ps -1 is needed for charmed hadrons of decay rates: T(A C ) = 
5 ps" 1 and T(D°) = 2.4 ps" 1 . 

The four-quark operators are estimated using phenomenological models. Hence their size is questionable. Never- 
theless, it is the only way now to do. The four-quark operators are related to the probability of finding the Qg pair 
at the origin simultaneously using quark models, denoted by 1^(0)1^. The wave function is evaluated relating it to 
the mass splitting of hadrons arising out of the heavy-light quarks interaction. In Ref . Q , Rosner evaluated the wave 
function, using the hyperfine splitting in a similarly heavy-flavoured baryon, utilising the DELPHI value on - 
splitting H [|] Neubert and Sachrajda [Q introduced hadronic parameters accounting for hybrid renormalisation while 
parameterizing the four-quark matrix for b-flavoured hadrons. The hadronic parameters are yet unknown. Their val- 
ues are obtained from QCD sum rules. As an extension to c-flavoured hadrons, Voloshin || studied the the charmed 
baryons. In this approach, the authors of Ref. || analysed the inclusive charmed-baryon decays to fix the hierarchy 
of charmed lifetimes. 

In the present study, we take into account the contribution coming from the exchange of hard gluons when two 
different scales, the heavy quark mass, m, and the QCD scale, fi, are involved, while estimating the wave function 
in quark model. We treat the coupling constant corresponding to a meson and a baryon differently. It is found that 
the wave function, and hence the enhanced decay rate, is large in both the cases of b and c and explains part of the 
lifetime differences between A& and B° as well as between A c and D°. 
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1 Though the DELPHI value has not hitherto been confirmed, its use does not alter the goal anyway since the central value of 
the mass splitting due to hyperfine interaction is expected to be around 50 MeV. 
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The enhancement in the decay rate of Af, is arising out of the processes involving four-quarks |MM-ttQl: (a) the 

weak scattering process bu >cd in the Af, involving matrix elements between hadronic states of (bb)(uu) operators 

and (b) the process contributing to the Pauli interference involving matrix elements of operators (bb)(dd) operators. 
Hence, the enhancement of the Af, decay rate is given by: 

Ar(A fc ) = ^\^(0)\l u \V ud \ 2 \V cb \ 2 m 2 b (l - x) 2 [C 2 _ - (1 + z)C + (C_ - C+/2)] (1) 

where x = m 2 /m 2 ; C_ and C+ = C 1 / 2 are the short distance QCD enhancement and suppression factors for quarks 
in a colour antitriplet and sextet respectively: 

C_ = [a s {ml)/a s {m 2 w )] A lf\ / 8 = ll-2n / /3 (2) 

where u/ is the active quark flavours between rah and raw- The C_ term corresponds to the weak sattering process 
bu — + cd — > bu and the other term represents the destructive interference between the two intermediate d quarks in 
the process bd — * cudd — ► cd. 

The wave function for the bu pair (in the initial baryon), |^(0)|L, in Eq. (1) is of the form: 

where 

AM(B ijk ) = ^o.^<^>|tl!r(0)|» (4) 
9 m,.m, J 

AM(M i3 ) . ^ a .<^>|.( ,| S (5) 

are the hyperfine mass splittings in a baryon and in a meson respectively. There is the colour factor 1/2 in the baryonic 
case due to the quark composition: a heavy quark and two light quarks. Under isospin symmetry, the effective masses 
of light quarks are equal. In the same token, wave functions for bu and bd pairs are equal. Equation (3) is obtained 
for the values of < Si.Sj > = (1/4, -3/4) with spin (0, 1) for the meson and < Si.Sj > = (1/4, -1/2) with spin (1/2, 
3/2) for the baryon with S qq = 1. In Eq. (3), £ is the ratio of the coupling constants governing a baryon and a meson. 
The coupling governing a baryon is stronger than that of a meson. Though £ has been chosen as a free parameter 
varying from 0.25 to 1.5, it can be exactly calculated JTlJ] . In rigorous sense, it should be more than unity. The 
wave function on the right hand side of Eq. (3) corresponds to the matrix element for B-meson decay into vacuum, 
parameterised as, 

< 0\q W Q\B > | 2 = f B M% (6) 

and to the relation obtained in the non-relativistic limit |l2j 

I < 0\<IWQ\B > | 2 = 12M B |v&(0)| 2 fl (7) 

Both Eqs. (6) and (7) characterise the same process involving a heavy quark and a light quark (and gluons) but are 
normalised at two different scales: the normalisation point of Eq. (6) is /i(«R _1 ) which is the order of the virtuality 
of the quarks inside the meson, and the scale of the Eq. (7) is the mass of the heavy quark. Therefore, in order to 
account for the hard gluon exchange between the light and heavy quarks, Eqs. (6) and (7) have to be related. The 
relation turns out to be an evolution equation of current operators p3jl4 j: 



and 



< 0|gVy 5 Q|fl > (mo) =< 0\q~f^ 5 Q\B > (^)[a s (^)/a s (m Q )]''/ /3 , fi < m Q (8) 

where 7 (= 2) is the hybrid anomalous dimension. The coupling constant defined through leading order renormalisation 
group equation can be expressed as 

a,(m) = ^ (9) 

I 1 - ^ a s(^)ln(fi/m)} 
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Using eq. (8), one can get the eq. (3) of the form 



|V " (0)L = 3 AM(M i j) 12 * I, 1 - iT^m J (10) 

The parameter £ describes the ratio of the coupling constants of baryon and meson. The other term in the bracket 
represents the logrithmic effects due to the hard gluon exchange that would be expected when one goes down to a 
scale as small as the hadronic scale from the heavy quark mass. This obviously modifies the expectation values of the 
wave function density at the origin. 

The choice of values of the parameters are given in Table I. The AM (Bgk) for A;, is given by M(E£)- M(5V) of 
DELPHI ||. For c, the AM(B ijk ) is given by M(E*) = 2517.5 1.4 MeV and M(E C ) = 2452.2 0.6 MeV flj. 

The enhanced decay rate is now given by Eq. (3) alongwith eq. (10). The results for the enhanced decay rate of 
Ab and A c are given in Tables II and III respectively. The wave function density is, roughly speaking, independent 
of the hadronic scale value. However, its dependence upon the parameter £ is important. For £ is graeter than one, 
the enhanced decay rate becomes closer to the enhancement required to explain the smaller lifetimes of Af, and A c 
baryons. 

In conclusion, it is demonstrated that the four-quark operators appearing at the third order in 1/m expansion 
explains the needed enhancement in the decay rates of Aq. Though subtle, the couplings of meson and Af, baryon 
make much difference. On the other hand, if the couplings are considered equal, the four-quark operators still account 
for the difference in lifetimes the A& baryon and B meson. 
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TABLE I. Choice of values of parameters El (n f = 4; Aqcd = 0.2 GeV). 



Q 


A U(B ijk ) 


A M(M i? ) 


M H 


fa 


mQ 


« s (m Q ) 


c+ 


C_ 


b 


0.056 


0.046 


5.2789 


0.190 


5.1 


0.20 


0.87 


1.31 


c 


0.0653 


0.142 


1.8693 


0.240 


4.7 


0.29 


0.79 


1.57 
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TABLE II. Enhanced Decay Rate for Af,. Needed Enhancement is 0.2 ps 





H = 1 GeV 




H = 0.5 GeV 






|*(0)| M - (10" 2 GeV 3 ) 


Ar(A f) )(ps- 1 ) 


l*(0)| M " (HT 2 GeV 3 ) 


Ar(A 6 )(ps- 1 ) 


0.50 


1.3154 


0.0763 


1.3422 


0.0779 


0.75 


1.9732 


0.1145 


2.0133 


0.1168 


1.00 


2.6309 


0.1527 


2.6844 


0.1558 


1.25 


3.2886 


0.1909 


3.3555 


0.1948 


1.50 


3.9463 


0.2291 


4.0266 


0.2337 



TABLE III. Enhanced Decay Rate for A c . Needed Enhancement is 2.6 ps 





H = 1 GeV 




At = 0.5 GeV 






|*(0)| cd - (10" 3 GeV 3 ) 


Ar(A c )(ps- 1 ) 


|\K(0)| cd - (10- 3 GeV 3 ) 


Ar(A c )(ps- 1 ) 


0.50 


2.7718 


1.2005 


2.8216 


1.2221 


0.75 


4.1578 


1.8008 


4.2325 


1.8332 


1.00 


5.5437 


2.4011 


5.6433 


2.4443 


1.25 


6.9297 


3.0015 


7.0542 


3.0553 


1.50 


8.3156 


3.6017 


8.4650 


3.6665 
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